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The dispersion properties of an anisotropic medium, such as the ionosphere in presence of the 
earth's magnetic field, are conveniently represented by the variations of the refractive index with the 
direction of propagation. It is shown that, for a number of cases of practical interest, the variation 
of the index n with the angle 6 takes a special form which can be expressed by a graphical construction; 
the result n is read on a circular scale, called the index circle, at the points where it is intersected by a 
straight line which depends in a simple manner on the direction of propagation 0. The same construc- 
tion applies when the wave is evanescent and then it gives the imaginary value of the refractive index 
or extinction index. 

The range of application of this construction includes the Appleton-Hartree equation, its generaliza- 
tion for any number of ions, the hydromagnetic waves in a medium with finite compressibility, and 
the various approximations corresponding to limiting or transitional cases — Astriim waves, uniaxial 
media, and Booker's hydromagnetic approximation. It could also be applied to magnetized ferrites. 

tor most applications the index circle representation can take the place of the dispersion surface. 
The main features of this surface — the number of sheets, whether they are real or imaginary, and the 
direction of the asymptotes — are all displayed on the index circle. If the need arises, it is easy to 
trace the dispersion surface and to follow its deformations as the parameters of the medium are varied. 



1. Introduction 

In analyzing problems of propagation, diffraction, 
or radiation in a uniform anisotropic medium it is 
sometimes helpful to consider first the plane waves 
that can exist in that medium. These waves, called 
characteristic waves, can then be added to form solu- 
tions of the wave equations that satisfy desired bound- 
ary conditions. The characteristic waves have specific 
polarizations and their wave vectors cannot be arbi- 
trary, but must have their extremities on a surface 
2, called the dispersion surface. This surface is 
also called the index surface because k/ko = n can be 
considered as the refractive index in the direction 
of vector k. Another surface, related to the index 
surface by inversion, gives the phase velocity c/n 
as a function of direction. 

The equation f{n, 6) = of the surface S can be 
found for various ranges of frequency and various 
conditions in the ionosphere. Even when this equa- 
tion is only of the second degree in n^ the dependence 
of the roots on the parameters of the medium, and on 
the angle 6, is complicated enough to make it difficult 
to visualize the variations of n with the parameters. 
Even to obtain a sketch of the dispersion surface re- 
quires considerable work. The purpose of this paper 
is to point out some similarities between various dis- 
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persion relations, in particular in the angular depend- 
ence of n, and to present a graphical construction 
which simplifies the evaluation of the refractive index. 
This construction was introduced by Deschamps and 
Weeks [1962] for a single component magnetoplasma. 
It is now extended to any gyrotropic medium such 
as a multicomponent cold magnetoplasma or a mag- 
netized ferrite. It is shown to apply also to the hydro- 
magnetic waves in a compressible fluid. 

2. General Form of the Dispersion Relation 

We shall consider exclusively dispersion relations 
of the form 

(L-L+)(L-L-)cos2^ + (L-Lo)(L-L')sin2 = 

(1) 

where L = n~^ and the coefficients L+, L-, Lo, L' are 
known functions of the frequency. The angle be- 
tween the direction of propagation and the z-axis is 0. 
The variable L, when it is positive, is proportional to 
the square of the phase velocity in direction 6. The 
coefficients L+ and L- are the values of L for ^ = 0, 
that is, for longitudinal propagation. The coefficients 
Lo and L' are the values of L for = 7t/2, that is, for 
transverse propagation. The reason for the pecuhar 
notation for the coefficients will become clear when 
considering a gyrotropic medium. The four numbers 
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L completely specify the surface S. Although the 
form of the relation (1) may seem very special, it does 
represent many cases of interest. Before proving this 
point by giving examples the construction which re- 
sults from (1) will be described. 

3. Description and Use of the Index-Circle 

The first step is to construct a circular scale on 
which the values of n will be read. This is done as 
shown in figure 1; a linear scale in L = n~^ is projected 
on the circle Y tangent to it using the bottom point 
of the circle as a center of projection. Alternatively, 
a linear scale in K = n^ can be projected from the top 
point, giving the same result. A graduation in n is 
shown in the figure, but a scale in phase velocity c/n 
could equally well be carried by the same circle. 
The points on the circle will be designated by the 
values of n, but sometimes by the corresponding val- 
ues of K or L. The points for two recriprocal values 
of n are on the same vertical. The bottom point of 
r is n = and the top point n = ^. Real values of n, 
corresponding to propagation, are represented to the 
right of the vertical axis. Imaginary values, corres- 
ponding to extinction,^ are represented to the left. 

Consider now the relation (1) and use ^ = sin^ ^ as a 
variable instead of 6. If the roots of (I), denoted by 
Li{t) and Luit) are plotted on the index circle, the line 
Dt joining the two representative points will depend 
on t. We note that Lj and Ln are in one-to-one cor- 
respondance because Lj defines t which defines L^. 
This correspondance is obviously algebraic and re- 
ciprocal hence it is an involution. From a well-known 
theorem in projective geometry it results that the line 
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Figure 1. Graduation of the index-circle. 

For illustration two projection lines are shown; from the top point to K='0A9 through 
n = 0.70 and from the bottom point to L = 0.25 through n = 2.00. In this figure the refrac- 
tive index is denoted by n = fi—jx while in the text fx denotes the permeability and x the 
susceptibihty. 



Figure 2. Variation of the line Dt with t = sin^ 0. 

The squares of n+, n , ru,, and nx are K+, K , Ko, and Ki; they are the reciprocals of L+, 
L , Lo, and L' . Because of the symmetry of the dispersion surface only one quarter of its 
meridional section is represented. The dotted line indicates an imaginary sheet of the 
surface corresponding to extinction. 



Dt goes through a fixed point J. We know two posi- 
tions of the hne Dt : Do which joins n+ to n- and Di 
which joins no to m. Thus the point / is known. Re- 
writing (1) in the form 

{L-L^){L-L-)-\-t(L^^L--Lo-L'){L-L2) = (2) 

shows that the roots for t = ^ are L = oo and L = L2 
= {L+L- — LoL')l(L+-\-L- — Lo — L')^ corresponding to 
n = and n = n2. These two points define Doc, which 
also goes through the point J as shown in figure 2. 
Finally, the dependence of Dt on t is one-to-one and 
algebraic. Thus the intersection of Dt with a parallel 
to Z)oo will describe a linear scale in i = sin^ 6. 

The procedure to find the refractive index in any 
direction is, therefore, as foUows. Mark on the index 
circle the values of the indexes for longitudinal and 
transverse propagation (fig. 2). Draw the lines Do 
and Di which intersect at some point /. Join this 
point to the bottom point n = to obtain D^o. Place 
a scale, linear in sin^ 6^ in such a way that it be parallel 
to Z)oc with its extremities and 1 respectively on Do 
and Di. Join the point J to the desired value of t on 
that scale. This gives the line Dt which intersects 
the circle F at Ui and Jiu. As the angle varies from 
to 7r/2, the line Dt sweeps through the shaded region 
in figure 2. In this instance one sheet of the disper- 
sion surface is real while the other is partly imaginary. 
The surface is a surface of revolution about the direc- 
tion of the constant magnetic field. A meridional 
cross section of this surface, deduced from the index- 
circle, is sketched in figure 2. The angle 6c of the 
asymptote is obtained by joining n = ^ to / and read- 
ing tc on the scale; here tc = 0A7, hence ^c = 43°. 
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4. Multicomponent Magnetoplasma 

The first application of relation (1) is to a cold mag- 
netoplasma with an arbitrary number of components. 
Each type of charged particle in the plasma is speci- 
fied by a number density ni, an electric charge qi with 
sign e/, and a mass m,. Two derived quantities are the 
plasma frequency [[i= {riiqjlniieo)^/" and the gyrofre- 
quency ftj^l/jioHoqilmil where Ho is the applied mag- 
netic field in the direction of the z-axis. Collisions 
are neglected and the susceptibility of the medium 
is the sum of the susceptibilities of the components 



x = Xx<" 



(3) 




The susceptibility matrix for the iih component of 
the plasma has eigenvectors e.s defined by 

e±= ^ (x=hyy), eo = z 

where x, y, and z are the reference unit vectors and 
the subscript s takes the values H-, — , and 0. The 
corresponding eigenvalues are 



n? 



X.V' 



(x){a) -\- S€jili) 



(4) 



The dielectric constant of the medium becomes a 
matrix K= 1+ X which has also the eg for eigenvectors 
and 



for eigenvalues. Thus K takes the special form 



(5) 



K 



K' jK" 
-jK" K' 
K„ 



(6) 



with K± = K'±K". This particular form of the di- 
electric matrix characterizes a gyrotropic medium. 
It is convenient to introduce the matrix L = K^' which 
has the same eigenvectors e,, as K with eigenvalues 
Ls^K~^. The matrix L takes the form 



L = 



L' 


jL" 





-JL 


L' 











U 



(7) 



with L+ =L' ±L". Let us denote by K\ the harmonic 
mean of K+ and K- and by L\ that of L+ and L-; ^i is 
the reciprocal of L' and L\ that of K' . 

With these notations, the dispersion relation takes 
precisely the form (1). However, while in (1) the four 
parameters L+, L_, Lo, L' were arbitrary, they must 
now be related since K depends only on the three 
numbers Ks-. The relation is that L' is the average of 
L+ and L-, hence K\ is the harmonic mean of K+ and 
K-. A construction of K\ and K' on the index circle 



Figure 3. Construction of the average K' and of the harmonic mean 
Ki o/K+ and¥>.-. 

They are obtained by projecting respectively the points n = <» and n — through the pole 
P of the Une K^K-. 



when K^ and K^ are given is shown in figure 3. Its 
justification is that both cross-ratios {K^K K'^) and 
(K+KKiO) are equal to — 1 and this relation is invar- 
iant under the projection that defines the circular 
scale. In figure 2 it has been assumed that K\ = nl 
is the harmonic mean of K+^n\ and K- = rii; thus it 
is an illustration valid for a gyrotropic medium. 

A magnetized ferrite is also a gyrotropic medium 
because its permeabihty matrix has the form (6). The 
dispersion relation takes also the form (1), but with 
^0=1. The resulting surfaces 2 will not be discussed 
here. 

5. Appleton-Hartree Dispersion Relation 

For frequencies much higher than the ion plasma 
frequency and the ion gyrofrequency, the properties 
of the ionosphere are well represented by the Apple- 
ton-Hartree equation. This equation assumes that 
the ions are at rest and takes into account only the 
motion of the electrons. The susceptibility X is that 
of the electrons. In this particular case there is one 
more relation between the three eigenvalues K.s, and 
it happens because Xo is the harmonic mean of x+ 
and X • This implies 



1 



+ - 



1 



Ko-l K^-l K--1 



(8) 



Thus K+ and K- are sufficient to determine Kq and Ki 
and therefore the whole dispersion surface. It is easy 
to verify that (8) implies that L2 = K2 = n2=l; the line 
Dec is now a fixed line joining the point n = io the 
point n=l. This particular relation was used by 
Deschamps and Weeks [1962], and the construction 
of Do and Di in terms of the magneto-ionic parameters 
X and Y was discussed in detail. 

Figure 4 shows a sketch of the index-circle represen- 
tation of the dispersion relation for various pairs of 
values of the parameters X=WIm^ and Y^ = fl^/(jt)^. 
This display is limited to a region in the vicinity of 
the plasma resonance and the gyroresonance of the 
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Figure 4. Index-circle representation of dispersion properties in a magnet op las ma. 

The vertical and horizontal axes are Y^^il^loj^ and X=\\^l(o^ respectively. The other lines in the (X—Y^) plane indicate where 
changes in the type of dispersion surface occur. The index no is distinguished by a small dot. 



electrons and to higher frequencies. The hnes in 
the (X—Y^) plane where the type of dispersion surface 
changes are indicated. This should be compared 
with plots given by Clemmow [1955] and Allis [1962] 
for the index surface and the phase velocity surface 
as a function of X and Y^. 

6. Hydromagnetic Approximation 

Consider a plasma made up of electrons and only 
one type of positive ions. The electrons are charac- 
terized by the frequencies lie and He, the ions by Hi 
and Hi. Because of the neutrality of the plasma, 



ni 



Hi 



(9) 



When the frequency co is much larger than Hi and Cli 
the Appleton-Hartree equation is justified. When 
o) is much smaller than He and He another simplifica- 
tion becomes possible [Booker, 1963]. This simpli- 
fication can be justified by finding approximations 
for the three eigenvalues Kg. Two of these, K+ and 
K-^ are 



K. = l- 



n^+n? 



(C0±110(c0=hae) 



(10) 



In the numerator Hi can be neglected compared to 
He and in the denominator o) can be neglected com- 
pared to fie hence, 



K^=^l-- 



m 



iieid) ± no 



using (9) 



K^^l + 



m 



[li{ni±a)) 



the third eigenvalue 



W H^ 



(11) 



(12) 



Thus the electron characteristics are eliminated and 
the subscript i will be dropped in what follows. 
Introducing 



KA = n\ = \^ 






(13) 



we see that the cross-ratio {K+K-Ka\) equals —1, 
hence the four corresponding points on the index- 
circle form an harmonic division. Consequently 
the tangents to F at point Ka and at point n = \ inter- 
sect on the line Do joining K^ to K- (fig. 5). The line 
Di joins Kq = ^ lo the point Ki constructed as in figure 
3. The lines Dq and Z)i intersect at 7, and D^ is the 
line from ai = to 7 . The variation of Dt with t = sin^ 
is obtained as in section 3 by drawing a parallel to 
/)oo and marking ^=0 and ^=1 on Dq and Di respectively. 

7. Astrom Waves 

If the frequency co is much lower than flj, a further 
simplification results: K+ = K- = Ka. The line Do 
becomes tangent to T at Ka and Di joins Ka to Ko=^ 
while Doo joins Ka to (fig. 6). The velocity correspond- 
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ing to ua reduces, when II ^ 11, t( 



c 

JlA 



'4-<T 



(14) 



where pi is the ion density rumi. This is the ionic 
Alfven velocity. 





Figure 5. Hydromagnetic approximation. 

The construction is for the case where llj and (1, have the same order of magnitude; tia 
is neither equal to 1 nor to <». The frequency w < 11; hence n~ is to the right of the vertical 
axis, otherwise it would he to the left of it. All points labeled n with some subscript could 
have been labeled K with the same subscript. 



Figure 7. Uniaxial medium. 

Three possible cases are shown with the corresponding dispersion surfaces: no real 
All, no real < n\, and ra© imaginary. The transition nfl = <» is the case of figure 6. 




Figure 6. Astrdm waves. 

Index-circle representation and dispersion surface when a> < O. 



The dispersion relation can be factorized. As Dt 
rotates about Ka one of the solutions is n = nA and the 
other is /i = n^/cos 6. This can be deduced from the 
general construction for Dt. The dispersion surface 
is made up of a sphere of radius Ha, corresponding to a 
wave with transverse polarization of the E vector, and 
of two planes z = ±nA, corresponding to a wave with 
transverse polarization of the H vector. 

8. Uniaxial Medium 

For a gyrotropic medium the dispersion relation is 
invariant under rotation about the z-axis and the matrix 
K has three distinct eigenvalues as explained above. 
When the two eigenvalues K+ and K- are equal, the 
matrix K becomes diagonal, ^" = 0, Ki = K\ and the 
medium is said to be uniaxial. The dispersion surface 
decomposes (fig. 7) into a sphere n = ni and a quadric of 
revolution with semi-axes rii in the z-direction and no 
in any direction of the xy plane. The Alfven-Astrom 
dispersion relation (fig. 6) is the special case where 
Aio = °o. 

9. Hydromagnetic Wave in a Compressible 
Fluid 

Let us consider a conducting fluid with finite com- 
pressibility in the presence of a constant magnetic 
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field Ho. The field quantities are the electric vector 
E, the magnetic vector H, and the velocity v. These 
are related by field equations as described in any work 
on magnetohydrodynamics. See, for example Jackson 
[1962] or Stix [1962]. Parameters describing the 
medium are the acoustic speed a= (ypo/po)^/^ and 
the Alfven speed A=Ho{iJiolpoy^^- In these formulas 
Po is the density, po the pressure, ixq the permeability, 
and y the ratio of specific heats. The formula for A 
has the same form as (14) with po, the fluid density, 
replacing the ion density p,. To the two character- 
istic speeds we can associate two refractive indexes 
Jia = cja and ra = cjA and the numbers La = n~^ and 

The dispersion relation is shown to decompose into 
two equations. The first one 



or 



L = La cos^O 



n = a/cos 6, 



(15) 



applies to the Alfven wave where the velocity v is 
transverse to the wave vector k and to the magnetic 
field Ho (shear mode). The corresponding dispersion 
surface is made up of the two planes z = ±nA. The 
second equation takes the special form (1) 



(L-La) (L-La) cos^d^L{L-LA-La) sin2^ = 



or 



(L-La) {L-LA)-LaLAsin^e = 0, 



(16) 



(17) 



The index-circle representation can be applied to this 
situation as shown in figure 8. Because the two pairs 
of roots La and La for ^ = 0, and and La-^LA for 
6 = 77 12 have the same average the lines Do and D\ 
intersect on the tangent to the circle F at the bottom 
point. This tangent is in fact D^o because (17) shows 
that the roots for s\n^d = ^ are both L=oo. Figure 8 
is for a> A and the resulting dispersion surface is 
drawn below the circle F. If a > ^ the surface would 
have the same shape but the planes z = ±nA would 
touch the open sheet instead of the closed one. 

When a is very small compared to A and to c, ria 
can be considered infinite and the compressibility can 
be neglected. The open sheet recedes to infinity and 
the closed sheet becomes a sphere. The situation is 
quite similar to that of the Astrom waves shown in 
figure 6. 

10. Conclusions 

The method described applies only to a special form 
of the dispersion relation. The medium must be loss- 
less and have symmetry of revolution. The equation 
must be, at most, of the second degree in n^ or decom- 
pose into factors having this property. These condi- 
tions, however, are met in a number of cases: magneto- 
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Figure 8. Hydromagnetic waves in a compressible fluid. 

The horizontal asymptote to the open sheet is given by /Cs = n^ = n^ + n^ and the finite 
value of the transverse refractive index by Lc = n^^ = n~^ + n~^. 



ionic theory of cold plasmas, hydromagnetic theory of 
compressible fluids, and theory of magnetized ferrites. 
The index-circle representation gives the imaginary 
values of n (extinction index) as well as the real ones 
for any direction of the wave normal. It shows at a 
glance the type of the dispersion surface; whether it 
has one or two real sheets, or none. It shows the zeros 
and infinities of the refractive index, cutoff, and reso- 
nances, and gives the critical angles for which they 
occur. 
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